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A new model of Laplacian stochastic growth is formulated 
using conformal mappings. The model describes two growth 
regimes, stable and turbulent, separated by a sharp phase 
transition. The first few Fourier components of the mapping 
define the web, an envelope of the cluster. The web is used 
to study the transition and the dynamics of large-scale fea- 
tures of the cluster characterized by evolution from macro- to 
micro-scales. Also, we derive scaling laws for the cluster size. 

PACS numbers: 64.60.Ak, 05.20.Dd, 41.10.Dq 



The relation of the two kinds of Laplacian dynamics, 
stochastic [Q| and deterministic ||^ , represents a very in- 
teresting problem. In the stochastic dynamics, such as 
diffusion- limited growth dielectric breakdown Q, or 
fracturing, the cluster grows by particles diffusing and 
sticking to it, or by emanating lightning strikes, or cracks, 
with probabilities determined from Laplace's equation. 
The growing cluster is a scale-invariant fractal object 
with non-trivial geometric characteristics that have been 
extensively explored by a combination of numerical and 
analytical methods • In deterministic growth , 

exemplified by the Hele-Shaw dynamics, the boundary 
of the growing object moves with a local velocity deter- 
mined by Laplace's equation. In this problem, attention 
was focused on the fingering instability, which at small 
capillary effects leads to fractal-like patterns resem- 
bling DLA cluster. To draw an explicit relation of the two 
large-scale structures, however, turns out to be difficult, 
even numerically. 

A natural way to proceed would be to study the con- 
tinuum limit of stochastic growth, which appears to be 
non-trivial. Taking DLA as an example [||, the arising 
difficulty is that the naive continuum limit, understood 
as taking the particle size to zero, gives Hele-Shaw dy- 
namics with zero surface tension [Q, leading to finite- 
time singularities. In that, Laplacian growth has a lot 
of similarity with turbulence. The finite-time singular- 
ities conjectured for Euler dynamics of an ideal fluid, 
are essential for understanding turbulent flow The 
large-scale spectrum of velocity fluctuations in the flow 
is determined by the instability cascade in inertial range, 
with viscosity and thermal noise being relevant only on 
microscale, where the singularities are resolved. Also, 
in parallel with energy conservation in the inertial range 
of Euler dynamics, it is found that all integrals of Hele- 
Shaw dynamics are conserved in DLA growth as well, 
up to small fluctuations arising on microscale p^ ] . From 
that, it is natural to conjecture that the large-scale DLA 
dynamics is equivalent to Laplacian contour dynamics. 



with the flnite-time singularities being resolved on mi- 
croscale due to noise. 

In order to explore this similarity, we propose a 
new class of models, continual rather than lattice-like, 
that provides an explicit connection between stochastic 
growth and contour dynamics. The growth is represented 
as a random sequence of conformal maps with memory, 
which facilitates numerical and analytical treatment. (It 
has been shown |^ that conformal maps simplify analysis 
of the mass distribution in the conventional DLA prob- 
lem.) In this model there are different parameter regimes 
analogous to the low and high Reynolds fluid dynamics: 
weakly stochastic, macroscopically stable growth, and 
noise-driven turbulent growth. The transition between 
the two regimes is sharp. By varying the "Reynolds" pa- 
rameter, one can model several known growth problems, 
such as dielectric breakdown and DLA. 
Stating the problem: One step of Laplacian growth in- 
volves attaching a new object to the cluster with a prob- 
ability determined by the solution of Laplace's equation 
V^u = 0, with M = at cluster boundary. The probabil- 
ity for the new object to appear within the interval dl of 
the boundary is given by dP = \Wu\dl. The size of the 
new object (quantifled by its area) can be 

a) constant, as in DLA; 

b) proportional to some power of local fleld, as in dielec- 
tric breakdown. For reasons made clear later we choose 
to write the power as a — 2, so that the area~ |Vm|"~^. 

Also, there is a freedom to choose different shapes 
for the object. To model breakdown and fracturing 
problems, where an individual growth step is a light- 
ning strike, or a crack, we consider growth that involves 
one-dimensional objects (hereafter called "strikes"). For 
diffusion-controlled growth, it is more natural to take new 
objects roughly equal in all dimensions (called "bumps"). 

In the canonical lattice model of dielectric breakdown 
[Q the size of a new object is flxed, but the sticking prob- 
abilities are given by a power law: dP = \\/u\^dl. To 
relate it with our model, one compares local area growth 
rates, and flnds rj = a — I. Given that, and assuming 
that universality classes are determined only by the local 
area growth law, we expect that the two models lead to 
the same macroscopic properties. Relation to the DLA, 
known to be a particular case of dielectric breakdown 
with T] — 1, is thus expected at a — 2. However, the 
relation of the lattice and the mapping models, although 
conflrmed by simulations, needs to be explored further. 
Mapping representation of growth: We use conformal 
mapping F{z) of the domain \z\ > 1 to the exterior of 
the growing cluster. (Uniqueness up to reparameteri- 
zation follows from the Ricmann mapping theorem.) To 



describe the change of the mapping F due to a single new 
object added to the cluster we use the mapping f\,e{z) 
that maps the domain 1^1 > 1 onto a sub-domain by at- 
taching a strike or a bump to the boundary at 2: = e** 
(see insets of Figs. , ). The strike-mapping fx.e{z) is 
e^^ f\{e~^^ z) where f\ is given by 
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where A is a parameter describing the size of the strike. 
Below, we use only mappings with A ^ 1, in which case 
the strike length is 2\/A + 0{X^/^). The bump- mapping 
can be chosen in several ways. We use 
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< a < 1. The growth is described by composition of 
mappings: 



(3) 



where the single step mappings /a^.s^ (z) can be or (||). 
Note that the order of the functions in the composition 
is reversed with respect to the growth time sequence. 

Now we determine the dependence of the parameters 
9i and Ai on the growth step. From confornial invari- 
ance of the two-dimensional Laplacian it follows that the 
random numbers 6i are uniformly distributed in the in- 
terval < 6 < 27r, and uncorrelated. For Xi, the power 
law relationship of the new object area and local field, 
|Vu|"-2 



-, gives 
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(To obtain (|) one just writes the Laplace's problem so- 
lution in terms of the mapping, u{z) = uoReln_F'(z).) 
Using (^ and the composition rule (^) one can grow the 
object without solving Laplace's equation at each step. 

Using (||) and (^), we simulate the growth. The exam- 
ples of growth patterns obtained for a = 0, 1.5, 2 are 
shown in Figs. , , . In the simulation, the data that rep- 
resent the structure are the set of angles 6i, and the set 
of Xi. For N growth steps, the algorithm complexity is 
iV^, which compares quite favorably to the canonical di- 
electric breakdown model (e.g., see and comes close 
to the best algorithms known for the DLA model. 

There is a special parameter value a = at which the 
model (H) , (^ is conformally invariant, which means that 
the stochastic dynamics commutes with arbitrary confor- 
mal transformations. As a result, at a = the growth 
does not have any "memory" : according to Eq. M) , all 
Xi are the same, and thus all mappings fxi^e^^z) are un- 
correlated. As a increases from to higher values, the 
memory effects become stronger and eventually lead to 
non-trivial dynamics. 

The web: To gain some insight, it is useful to consider 
power series expansions in z~^ for Fn{z) and f\„,s,^{z): 
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asymptotically equal to z + Xnz{z + e ")/(z 
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the expansion for / we keep only low order terms in A„, 
since we are going to work with small A„. Although to 
linear order the cutoff in (^ is infinite, the finite cutoff is 
written in explicitly to indicate that / has some smallest 
scale, beyond which it is smooth (see below). 

By truncating the series ^ one obtains a function that 
maps the domain \ z\ > 1 to an envelope of the grown clus- 
ter (we call it the web). The web, even for a relatively 
small number of terms, is an accurate representation of 
the cluster (see Fig. ). The web is sensitive to the struc- 
ture of the outer growing part of the cluster, rather than 
to the stationary inner region, which makes the web a 
useful tool. One observes a change in the web grown at 
different a. For < a < 1 the web is rounded, becom- 
ing more circular at large times. On the other hand, at 
1 < a < 2 the web is rough at all times, the roughness 
scaling as the cluster size. 

An alternative way of imaging the growing cluster used 
in the literature is based on the averaged occupancy [Q. 
Its relation to the web represents an interesting problem. 
The web dynamics: Next, we define web dynamics. To 
obtain the web for a given cluster requires a numerical 
Fourier transform. Using the expansions (|), (§) and the 
rules (IH), (^), we can instead directly evolve the power 
series, calculating new coefficents of the series from old 
coefficients. To do this, (||) must be cutoff at some fixed 
number of terms. As long as this cutoff is greater than 
the cutoff in (^), the web dynamics yields good results. 
Even for a relatively small number of terms (say, 20 to 
50), it agrees acceptably with the exact dynamics, both 
in the size of the cluster taken as function of time, and 
in the overall shape. Several examples of grown webs are 
shown in the insets of Fig. . Details on the web dynamics 
algorithm will be published elsewhere ||l^ . 

One advantage of the web dynamics is that it takes 
much fewer computer steps than exact dynamics. Also, 
the web can be used to study the large-scale limit of the 
growth dynamics. If for any reason one is only interested 
in macroscopic characteristics, such as the overall growth 
rate, or the roughness of the cluster boundary, then the 
web is a more appropriate instrument than exact map- 
ping. Finally, the web provides a bridge between exact 
dynamics and its continuum limit (see below). 
Taking continuum limit. Consider a continuum limit, in 
which the attached objects are infinitesimally small. 
Then, the growth is described by deterministic equa- 
tion written for F{z). To derive it we can use the 
power series (||). First, substitute the expansion (^) of 
f\„.en(z) in the recursion relation (^), and expand in A„: 
Fn+i{z) = F„(z) + SFniz), where 



dF„[z) = — A„z- 



dz 
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Then the continuum hmit is implemented by averaging 
over 9n- Taking the integral over ddn and substituting 
An = AoI-FjI"") gives Shraiman-Bensimon equation 

F{z) = XozFM <f |i.,(e'«")|-i±f;!:^ (8) 
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For a — 2, which corresponds to constant area of at- 
tached object, i.e., to DLA, we recover the Hele-Shaw 
problem which has rich analytic properties pO[ | and can 
be solved in terms of poles of F(z) moving within the disk 
I z I < 1 . It turns out that generic dynamics leads to singu- 
larities occurring at finite time: F{z,t) ^ A{t)/{z-~w{t)), 
\wit to)\ ^ 1. 

Let us compare this behaviour with the web dynamics. 
The singularities of the dynamics (H) can be recognized 
in the cusp-like features appearing on the web (see Fig. ). 
After appearing, they become sharper, but later, instead 
of developing into real singularities, they just disappear 
being gradually replaced by other similar features nearby. 
It is easy to draw a parallel with the currently accepted 
picture of turbulence jll]] . In the inertial range, the fluid 
dynamics gives rise to singular vortex filaments. In the 
absence of viscosity these would develop into singular- 
ities. However, on the viscous range the singularity is 
suppressed, and thus a turbulent flow exhibits a random 
sequence of singular-like structures replacing each other. 

In the web, the analog of the spatial scale is the power 
of z. It can thus be said that sharpening of the cusps on 
a web corresponds to the shift to smaller scale. Cusps are 
suppressed at microscale set by the largest power of z in 
(^. This shortest scale is determined by the value of A, 
which corresponds to the amount of noise; in the contin- 
uum limit the noise vanishes and the microscale becomes 
infinitely small. In the case of web dynamics, the trun- 
cation of (H) may also fix a microscale. In both cases, for 
the exact mapping and for the web, we have singular-like 
structures developing in the "inertial range" , where the 
continuum equation (^) holds, evolving to shorter scales, 
and disappearing there. It is thus appropriate to call the 
long-time dynamics turbulent. 

The role of the noise due to the particle discreteness is 
unusual. The dynamics without noise would by itself pro- 
duce chaos, and even singularities. The noise suppresses 
singularities, and makes the system run forever, similar to 
the role of viscosity in a turbulent flow. Such behaviour 
can be contrasted to other stochastic growth problems, 
such as the Kardar-Parisi-Zhang dynamics |14], or ki- 
netic roughening of stable Laplacian fronts |15], where 
all stochastic properties arise due to the noise. In such 
problems, naive continuum limit gives rise to stable dy- 
namics, and thus is non-problematic. 
Phase transition. As function of a, there is a transition 
at a = 1 from stable to turbulent growth. The stable 
growth at a < 1 is similar for the strike and the bump 



models: the degree of roughness of the cluster boundary, 
scaled in the cluster size, decreases with time. 

The transition at a = 1 is sharp. To verify this, we 
plot the mean square of the fluctuations of the boundary 
rescaled by the object size (see Fig. ). Here, we grew 
the object on a periodic strip of finite width instead of 
the circular geometry discussed above, and used the first 
Fourier component of the mapping as a measure of fluc- 
tuations. As the bump size Aq becomes smaller, the fluc- 
tuations decrease at a < 1, indicating convergence to the 
continuum limit (^, but remain finite at a > 1. 

These results may not exactly apply to the circular 
cluster growth, due to certain differences in DLA growth 
in the two geometries. In contrast with the circular 
growth, in the long time limit, the growth in the strip 
geometry is statistically unchanging, which is why this 
geometry was selected for the phase transition study. 

The turbulent growth at a > 1 is sensitive to the 
model. In the strike-models the object grows by emanat- 
ing long chains of strikes containing almost all the den- 
sity. On the other hand, for the bump-models, at a > 1 
the growing object forms a self-similar fractal cluster, 
with the fractal dimension varying with a (see Figs. , ). 
The self-similarity was checked by using gyration radius, 
which is found to be a power law function of time over 
about 10^ steps, and by calculating fractal dimension us- 
ing the box counting method. The box counting fractal 
dimension was slightly lower, but was consistent with 
the gyration radius power law. (A similar discrepancy 
has been noted before The cluster grown at a = 2 

has properties identical to that of a DLA cluster. This 
is expected, since at a — 2 the new object size remains 
roughly constant throughout the growth. So, the a = 2 
bump models are equivalent to DLA. 
Scaling laws. Given the self-similarity of DLA, it is of 
interest to study scaling of the power series coefficients. 



For that we use the expansion of F^ 
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the analytic function whose absolute value gives the local 
electric field. The mean square of g„'s are plotted against 
n on a log- log plot (see inset to Fig. ). A renormaliza- 
tion group theory developed elsewhere [T9| ] predicts the 
slope of —2/5, also shown in the plot. The higher Fourier 
coefficients have a large effect on the growth rate of the 
cluster, as discussed below. 

The cluster size given by a^"^ in (||), by using (^), can 

be written in closed form: a^"'' = J|(l -I- Xj), < j < 
n. For a < 1, the object is roughly circular, and Xj 
is determined by oi. Noise in the higher components 
leads to some renormalization of Aq Aq , and then A^ = 



Xoia\ 



This equation has a power law solution: 
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The same power law exponent can be found from the 
circular-symmetric solution to Eq. (|^). 

In our simulation, at any a the growth of oi is indeed 
described by a power law, and at a < 1, the power law 



exponent is equal to 1/a, as predicted by Eq. (j^). How- 
ever, at a > 1, the exponent is somewhat larger than 1/a, 
and the deviation increases at larger a. This correction 
is due to the growth of the higher Fourier components. 
Any given one of the higher Fourier components eventu- 
ally stabilizes at some average value, but as the object 
grows, higher and higher Fourier components contribute 
to A, and the renormalization of Aq in Eq. (j^) increases 
with time, leading to a correction in the exponent. 

Even at large a it is still true that, on average, 
dai/dt = aiA, and so, since oi follows a power law in 
time, the average of A scales as i"^. This is equivalent to 
the electrostatic scaling law derived by Halsey (l6| 

In summary, we studied analytic models of Lapla- 
cian stochastic growth, formulated in terms of conformal 
mappings, and characterized two growth regimes: sta- 
ble (a < 1) and turbulent {a > 1), the transition at 
a = 1 being sharp. By truncating the series expansion 
of the mapping, we defined the web, geometric envelope 
of growing cluster. The web dynamics is used to demon- 
strate a relation between transient features of the turbu- 
lent growth and finite-time singularities of the determin- 
istic problem. The turbulent growth bump-models are 
found to be equivalent to the dielectric breakdown and 
DLA models. 
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Cluster grown using strike mappings with a = 0. Inset: 
A strike map applied to three circles of varying radii. 



Cluster grown using bump mappings with a = 1.5. Inset: 
A bump map applied to three circles of varying radii. 



Cluster grown using bump mappings with a — 2, and its 
web obtained by truncating the cluster mapping series 
expansion at 40 terms; Insets 1, 2 : Examples of growth 
using the web dynamics with 20 and 40 terms; Insets 
3, 4 ■ For comparison, the webs of the clusters grown 
using exact mapping obtained by truncating at 20 and 
40 terms. 
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Boundary roughness measured by the mean square value 
of the first Fourier coefficient for three different values 
of Aq. Inset: Mean squares of Fourier coefficients of F^^ 
scaled by the cluster size ai, averaged over 50 runs. The- 
oretical line with the slope —2/5 is drawn. 
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